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We report a search for the rare decays B'^ — > D^ K" and -B^ -^ D^ K*" in an event sample of 
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INTRODUCTION 



Charged B meson decays in which neither constituent 
quark appears in the final state, such as B'^ — >■ D'^K^*'^, 
are expected to be dominated by weak annihilation di- 
agrams with the bu pair annihilating into a W'^ boson. 
Such processes therefore can provide insight into the in- 
ternal dynamics of B mesons, in particular the overlap 
between the b and the u quark wave functions. Annihila- 
tion amplitudes cannot be evaluated with the commonly- 
used factorization approach [l|. As a consequence, there 
are no reliable estimates for the corresponding decay 



rates. Annihilation amplitudes are expected to be pro- 
portional to /s/ms where mB is the mass of the B me- 
son and fs is the pseudoscalar B meson decay constant. 
The quantity fs represents the probability amplitude for 
the two quark wave functions to overlap. Numerically, 
fB/iTT-B is approximately equal to A^, where A is the sine 
of the Cabibbo angle [ll y]- In addition, these ampli- 
tudes are also suppressed by the CKM factor | Vub | ^ A'^ . 
So far, there has been no observation of a hadronic B 
meson decay that proceeds purely through weak annihi- 
lation diagrams, although evidence for the leptonic decay 
B ^> TV has been found [3|. In theoretical calculations 



of nonleptonic decays, the assumption is often made that 
these amphtudes may be neglected. 

Some studies indicate that the branching fractions of 
weak-annihilation processes could be enhanced by so- 
called rescattering effects, in which long-range strong in- 
teractions between B decay products, rather than the 
decay amplitudes, lead to the final state of interest [3- 
Figure [T] shows the Feynman diagram for the decays 
B+^D+K'-*'>° and B+^D+Tr° 0], and the hadron-levcl 
diagram for the rescattering of Dfn'^ into D'^K'-*^^. Sig- 
nificant rescattering could thus mimic a large weak anni- 
hilation amplitude. It has been argued [2] that rescat- 
tering effects might be suppressed by only A"*, com- 
pared to A^ for the weak annihilation amplitudes, ren- 
dering the 5+ -^ D^K^*'>^ decay rate due to rescat- 
tering comparable to the isospin-related color-suppressed 
_B" -^ £)0/i'(*)o decay rate of approximately 5 x 10~^. 



resonance is used to study backgrounds from e+e~ -^ qq 
{q = u, d, s, or c) processes, which wc refer to as contin- 
uum events. 

The BaBAR^ detector is described in detail elsewhere 
[lO| . Charged-particle tracking is provided by a five layer 
silicon vertex tracker (SVT) and a 40 layer drift cham- 
ber (DCH). In addition to providing precise position in- 
formation for tracking, the SVT and DCH measure the 
specific ionization, which is used for particle identifica- 
tion of low-momentum charged particles. At higher mo- 
menta {p > 0.7GeV/c) pions and kaons are identified 
by Cherenkov radiation detected in a ring-imaging de- 
vice (DIRC). The position and energy of photons are 
measured with an electromagnetic calorimeter (EMC) 
consisting of 6580 thallium-doped Csl crystals. These 
systems are mounted inside a 1.5 T solenoidal super- 
conducting magnet. Muons are identified by the instru- 
mented magnetic-flux return, which is located outside the 
magnet. 
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FIG. 1: Annihilation diagram for the decay B+^D+ZC'*'" 
(top). Tree diagram (bottom left) for the decay B'^ — > Dfiv'^ 
and hadron-level diagram (bottom right) for the rescattering 
contribution to B+^D+K'-'''^ via B+-^Dfn°. 



B'^ — > Z?+i^(*)° decays are also of interest because 
their decay rates can be used to constrain the an- 
nihilation amplitudes in phenomenological fits [l|, |5[. 
This allows the translation of the measurements of the 
5+— >-Z?°if ^*)"'' amplitudes into estimations of the \Vub\ 
suppressed amplitudes B^^D^K^*^^ %%■ None of the 
modes studied here has been observed so far, and a 90% 
confidence level upper limit on the branching fraction 
B{B+^D+K°) < 5 X 10"^ has been established by 
BaBAR 0. No study of B+ -^ D+K*° has previously 
been published. 

The results presented here are obtained with 426 fb^^ 
of data collected at the T{4S) resonance with the BaBAR 
detector at the PEP-II asymmetric e+e~ collider [3| cor- 
responding to 465 X 10^ BB pairs (Ng^). An addi- 
tional 44.4 fb~ of data ( "off- resonance" ) collected at 
a center-of-mass (CM) energy 40 McV below the T(4S') 



EVENT RECONSTRUCTION AND SELECTION 

The event selection criteria are determined using 
Monte Carlo (MC) simulations of e+e" -^ T{4S) -^ BB 
("BB" in the following) and continuum events, and the 
off-resonance data. The selection criteria are optimized 
by maximizing the quantity S/y/S + B, where S and 
B are the expected numbers of signal and background 
events, respectively. We assume the signal branching 
fraction to be 5 x 10~^ in the optimization procedure. 

The charged particle candidates are required to have 
transverse momenta above 100 MeV/c and at least twelve 
hits in the DCH. 

Candidate D+ mesons are reconstructed in the D+ — > 
K^Tr+n+ {Kinr in the following), D+ -^ A'!?7r+ (AT^tt), 
D+ -^ K-TT+TT+TT^ {Kmrn^) and D+ -)■ K°tt+7t° 
(ATgTrTT^) modes for the 5+ — > D'^K^ decay channel 
(DK). Only the first two modes are used for the B^ ^■ 
]j+X*^ decay channel {DK* in the following) since we 
find that including the Kmnr'^ and K'^tttt'^ modes in this 
channel does not appreciably improve the sensitivity of 
the analysis. 

The D~^ candidates are reconstructed by combining 
kaons (either charged or neutral depending on the chan- 
nel) and the appropriate number of pions. The charged 
kaons used to reconstruct the D'^ and A'*" candidates are 
required to satisfy kaon identification criteria obtained 
using a likelihood technique based on the opening angle 
of the Cherenkov light measured in the DIRC and the 
ionization energy loss measured in the SVT and DCH. 
These criteria are typically 85% efficient, depending on 
the momentum and polar angle, with misidentification 
rates at the 2% level. Kaons and pions from D decays 
arc required to have momenta in the laboratory frame 
greater than 200 McV/c and 150 MeV/c, respectively. 



The reconstructed D"*" candidates are required to satisfy 
the invariant mass (Md) selection criteria given in Ta- 
ble D 

The K^ candidates are reconstructed from pairs of 
oppositely-charged pions with invariant mass within 5- 
7 MeV/c^ of the nominal K^ mass ll|. This mass cut 
corresponds to 2-2.8 standard deviations of the experi- 
mental resolution and varies slightly among channels due 
to the different amounts of background per channel. For 
the prompt K^ candidates from the B'^ — >■ D^K^ decay, 
we require 1 — cosQ;^o(i3+) < 10~^, where aK°{B^) is 
the angle between the momentum vector of the K^ can- 
didate and the vector connecting the B^ and K^ decay 
vertices. For K^ daughters of a D^ decay, we require 
1 — cosa^o (-D"*") < 10~^, where a^o [D^) is defined in a 
similar way. 

The 7r° candidates are reconstructed from pairs of pho- 
ton candidates each with an energy greater than 70 MeV, 
and a lateral shower profile in the EMC consistent with 
a single electromagnetic deposit. These pairs must have 
a total energy greater than 200 MeV, a CM momentum 
greater than 400 MeV/c, and an invariant mass within 
lOMeV/c^ (for the KinnT^ mode) or 12MeV/c2 (for the 
K^tttt'^ mode) of the nominal 7r° mass [ll| . 

The K*'^ candidates arc reconstructed in the decay 
channel K*'^ — )- K'^'k~ . These charged tracks are con- 
strained to originate from a common vertex. The recon- 
structed invariant mass, whose width is dominated by the 
K*^ natural width, is required to lie within 40 MeV/c^ of 
the nominal K*'^ mass Uj. We define O-a as the angle be- 
tween the direction of flight of the charged K and the di- 
rection of flight of the B in the K*^ rest frame. The prob- 
ability distribution of cos 0h is proportional to cos^ 6'h for 
longitudinally polarized K*^ mesons from B — ?> DK*'^ 
decays, due to angular momentum conservation, and 
is approximately flat for fake (random combinations of 
tracks) or unpolarized background if*" candidates. To 
suppress fake and background K*^ candidates we require 
I cos 6'h I > 0.5. 

The B'^ candidates are reconstructed by combining 
one D+ and one K^ or K*'^ candidate, constraining them 
to originate from a common vertex. The probability dis- 
tribution of the cosine of the B polar angle with respect 
to the beam axis in the CM frame, cos 9b, is expected 
to be proportional to 1 — cos^ 9b ■ Selection criteria on 
I cos 9b I are channel dependent and are summarized in 
Table H 

We measure two almost independent kinematic vari- 
ables: the bcam-cncrgy substituted mass toes = 
y/{E^'^/2+po -phY/El -pb'^, and the energy differ- 
ence AE = E^ — Eq/2, where E and p are energy and 
momentum, the subscripts B and refer to the candidate 
B and to the e+e" system, respectively, and the aster- 
isk denotes a calculation made in the CM frame. Signal 
events are expected to peak at the B meson mass for toes 
and at zero for /S.E. Channel-dependent selection crite- 



ria on |A£'| are given in Table HI We retain candidates 
with Toes in the range [5.20, 5.29] GeV/c^ for subsequent 
analysis. 

In less than 1% of the cases, multiple J5+ candidates 
are present in the same event, and in those cases we 
choose the one with the reconstructed D^ mass closest to 
the nominal mass value Uj . If more than one B'^ candi- 



date shares the same D'^ candidate, then we choose the 
B^ candidate with AE closest to zero. 



BACKGROUND CHARACTERIZATION 

After applying the selection criteria described above, 
the remaining background is composed of non-signal BB 
events and continuum events, the latter being the dom- 
inant contribution. Continuum background events, in 
contrast to BB events, are characterized by a jet-like 
shape, which can be used in a Fisher discriminant J^ [12| 
to reduce this background component. The discriminant 
J^ is a linear combination of four variables trained to 
peak at 1 for signal and at —1 for continuum background. 
The first variable is the cosine of the angle between the 
B thrust axis and the thrust axis of all the other re- 
constructed charged tracks and neutral energy deposits 
(rest of the event), where the thrust axis is defined as 
the direction that maximizes the sum of the longitudi- 
nal momenta of all the particles. The second and third 
variables are the event shape moments Lq = "^iPi, and 
L2 = X^iPil cos^ip, where the index i runs over all tracks 
and energy deposits in the rest of the event; pi is the mo- 
mentum and 9i is the angle with respect to the thrust 
axis of the B candidate. These three variables are calcu- 
lated in the CM. Finally we use |Ai|, the absolute value 
of the measured proper time interval between the two 
B decays [8|. It is calculated using the measured sepa- 
ration along the beam direction Az between the decay 
points of the reconstructed B and the other B and the 
Lorcntz boost between the laboratory and CM frames. 
The other B decay point is obtained from the tracks 
that do not belong to the reconstructed i?, with con- 
straints from the reconstructed B momentum and the 
beam-spot location. The coefficients of F, chosen to max- 
imize the separation between signal and continuum back- 
ground, are determined with samples of simulated signal 
and continuum events, and validated using off-resonance 
data. We denote two regions: the fit region, defined as 
5.20 < Toes < 5.29GeV/c2 and -5 < J" < 5, and the 
signal region, defined as 5.27 < toes < 5.29 GeV/c^ and 
0< J"< 5. 

To reduce the importance of the continuum back- 
ground in the final sample we divide the events according 
to their flavor-tagging category [8| . We define the follow- 
ing exclusive tagging categories: 

• lepton category, events contain at least one Icpton 



TABLE I: Main selection criteria used to distinguish between signal and background events. Md,pdg is the nominal mass of 
the D+ meson [Tl. 



Selection criteria 



B^ 



"D^W 



B^ 



D+K* 



TT 



KlVTV 



Ktt-ktv 



XT 



ro- 



Ktvtt 



K^ 



|Mo,pdg| (MeV/c^) <12 (~ 1.9o-) <18 (~ 1.5a) <14 (~ 1.6cr) <22 (~ 1.6a) 
cosbsl <0.76 <0.77 <0.87 <0.85 

\/\E\ (MoV) <20 (~ 1.3a) <23 (~ 1.5a) <25 (~ 1.5a) <24 (~ 1.5a) 



<10 (~ 1.6a) <10 (~ 1.8a) 

<0.82 <0.84 

<19 (~ 1.3a) <19 MeV (~ 1.3a) 



TABLE II: Reconstruction efficiencies and expected numbers of events in the fit and signal region assuming B{B^ — >■ D'^ K^) 
B{B+ -^ D+K*°) = 5 X 10"®. 





region ■ 




B+ 


^D+7^" 




B+^ 


D+A"" 




Ktttt 


K-K-rni" 


Ki^ 


K'W 


Knn 


^STT 


Signal efficiency 


fit 
signal 


18.4% 
12.4% 


5.2% 
3.8% 


21.3% 
14.7% 


6.2% 
4.9% 


10.6% 
7.6% 


10.5% 
7.4% 


Signal 


fit 
signal 


14.1 ± 0.2 
9.6 ± 0.2 


2.5 ± 0.1 
1.8 ± 0.1 


1.81 ± 0.03 
1.21 ± 0.03 


2.4 ± 0.1 
1.9 ± 0.1 


15.8 ± 0.3 
11.3 ± 0.3 


1.70 ± 0.04 
1.20 ± 0.03 


Combinatorial BB background 


fit 
signal 


67 ±4 

7± 2 


157 ±4 
20 ± 2 


12 ± 2 
3 ± 1 


36 ± 3 

8 ± 2 


400 ± 10 
30 ±2 


42.8 ± 4 
6.4 ± 1 


Peaking BB background 


fit 
signal 


2.0 ± 0.2 
0.3 ± 0.1 


3.3 ± 0.4 

1.0 ± 0.2 


1.1 ± 0.2 
0.3 ± 0.1 


1.8 ± 0.5 
0.6 ± 0.2 


26 ± 2 
5.4 ± 1 


2.4 ± 0.3 
0.7 ± 0.2 


Continuum background 


fit 
signal 


2840 ± 40 
63 ± 6 


4860 ± 50 
104 ± 8 


640 ± 20 
12 ± 3 


1600 ± 30 

45 ± 5 


6100 ± 50 
129 ± 8 


630 ± 20 
13 ± 3 



in the decay of the other B meson; 

• kaon category, events contain at least one kaon in 
the decay of the other B meson, which do not be- 
long to the first category; 

• other category contains all the events not included 
in the two previous categories. 

The first two categories are expected to be less contam- 
inated by continuum background. We fit all three cate- 
gories simultaneously. Studies of simulated events show 
that using the tagging categories reduces the statisti- 
cal uncertainty on the measured branching fraction for 
the Ktttt mode by 5%, but leads to little gain for the 
other modes (which are less statistically significant them- 
selves). Hence, we use tagging information only for the 
Kinr channel. 

The BB background is divided into two components: 
non-peaking (combinatorial) and peaking. The latter can 
occur when one or several particles of a background chan- 
nel are replaced by a low momentum charged tt^ and the 
resulting candidate still contributes to the signal region. 
The largest contributions to the BB peaking background 
for the B'^ -^ D^K^ channel arise from the following 
decays: B^ -^ D^p^ with D+ decaying into signal chan- 
nels, B" -^ D"K" and B" -> D*°K°. To further reduce 
the contribution from the i?° — >■ D~^p~ background, the 
variable | cos^^o | has been introduced, where Oj^o is the 
K'^ helicity angle, i.e., the angle between one of the two 
pious from the K'^ and the £)+ in the K'^ rest frame. We 



reject events with | cos 9j^o \ greater than 0.8 for the Ktttt 
mode and 0.9 for all other modes. Based on MC studies, 
we expect no more than 1 BB peaking background event 
per mode in the signal region, after applying all selec- 
tion criteria (see Table HI)) . A similar study is performed 
for the 5+ — > D^K*^ decay modes. The main peaking 
backgrounds arise from B° -^ D+p^, B^ -^ D^K*~, 
and 5° — )- D"'"af. In all cases, the D^ decays into the 
signal decay modes. The number of BB peaking back- 
ground events expected in the signal region for the DK* 
mode are shown in Table |lTl 



Charmless B decays may also contribute to the peak- 
ing background. These decays can produce tt and K 
mesons with characteristics similar to those of signal 
events without forming a real D meson. The charmless 
background is evaluated from data using the D^ side- 
bands: events are required to satisfy the criteria 1.774 < 
Md < 1.840 GeV/c^ or 1.900 < Md < 1.954 GeV/c^. We 
obtain —1.7 ± 1.0 events for DK decays and —0.7 ± 2.1 
events for DK* decays. We estimate the charmless peak- 
ing background contribution to be negligible and assign 
a systematic uncertainty based on this assumption. 



The overall reconstruction and selection efficiencies for 
signal events, as well as the numbers of expected events 
for each background category, are given in Table [ill 



TABLE III: Expected errors on the branching fractions from toy MC studies depending on the branching fractions generated. 
The combined errors are obtained as results of hkehhood combination per each toy (see text for details). All the numbers are 
given in units of 10~*. 



Decay mode 


B 


= 5 


e=o 


Mean error [95% range] 


Mean error [95% range] 


B+ -^ D+K" 


Knn 


+3.3 
-3.0 


[2.7,4.0] 
[2.2,3.6] 


+2.8 
-2.4 


2.2,3.6] 
1.6,3.2] 


Knmr'-' 


+20 

-17 


[14,25] 
[10,23] 


+ 19 
-17 


13,24] 
9.4,22] 


K^n 


+12 
-8 


[7.3, 16] 
[4.6, 14] 


+11 
-8 


7.1,16] 
4.5, 14] 


RO-kttO 


+14 
-12 


[8.9, 18] 
[6.2, 16] 


+13 
-11 


8.3, 17] 
5.6, 15] 


combined 


±2.9 


[2.1,3.6] 


±2.5 


1.5,3.2] 


B+ -> D+K*" 


Knn 


+3.5 
-3.2 


[2.5,4.0] 
[1.8,3.6] 


+3.3 

-2.8 


[2.5,4.0] 
[1.6,3.8] 


Kin 


+15 
-11 


[9.8, 19] 
[5.8,16] 


+14 
-7.7 


[7.9, 17] 
[3.8, 14] 


combined 


±3.3 


[2.1,4.2] 


±3.0 


[1.8,3.9] 



TABLE IV: Branching fraction fit results in units of 10 ^, with statistical uncertainties. A'^^ are the yields of the fitted species, 
and B represents the calculated branching fraction for each channel. 



Decay mode 



A^si 



N 



BB 



Nc, 



B 



B+ -i- D+K° 



Knn 
Knni\ 
Kin 
Kinn' 
combined 







-11-9 + ^:1 

10 

9 

5.3 

4.5 

4.5 



10 
0.6 



-6.7 ; 



2.8 



70 ± 27 2690 ± 57 

111 ± 51 6516 ± 94 

20 ± 14 381 ± 23 

36 ± 22 1270 ± 41 



+ 2.4 

- 2.0 
+ 20 

- 17 
+ 15 

- 13 
14 + 9.2 



-4.2 
20 
0.7 



-3.4 



6.2 

+ 2.2 
- 1.8 



B+ ->■ D+K*° 


Knn 
Kin 
combined 


-15.6 1«:I 

-11.4 t 3.5 


463 ± 63 6338 ± 98 

35 ± 15 547 ± 27 


5 + 2.9 
ti.U _ 2.1 

-33 t ^°o' 

= 0+2.3 
-O.d _ 2.0 



FIT PROCEDURE 

The signal and background yields are extracted by 
maximizing the unbinned extended likelihood 



£ = (e-^7^!) 



N 



■ N'''-l[f{^,\e,N'). 



(1) 



Here Xj = {tties;-^}, ^ is a set of parameters, N is the 
number of events in the selected sample, N' is the expec- 
tation value for the total number of events, and 



fi^\0,N')^ 



iVsig/.ig(x|0)+E.iVB./B.(x|g) 

N' 



(2) 



with /sig(x|6') and /^^(xl^) the probability density func- 
tions (PDFs) for the hypothesis that the event is a signal 
or a background event, respectively. The Bi are the dif- 
ferent background categories used in the fit: continuum 
background; combinatorial BB background, and peaking 



BB background. A^sig is the number of signal events, and 
Nsi is the number of events for each background species 
B,.' 

The individual probability density functions are de- 
fined by the product of the one-dimensional distributions 
of m-ES and J^. Absence of the correlations between these 
distributions is checked using the MC samples. The sig- 
nal mES distribution is modeled with a Gaussian func- 
tion. The continuum and non-peaking BB background 
Toes distributions are modeled with two different thresh- 
old ARGUS functions defined 131 as follows: 



A{x) 




(3) 



where xq represents the maximum allowed value for the 
variable x and c accounts for the shape of the distribu- 
tion. The Toes distribution of the peaking BB back- 
ground is modeled with a Grystal Ball (GB) function 



14| . The CB function is a Gaussian modified to include a 
power-law tail on the low side of the peak. The T distri- 
butions are modeled as the sum of two asymmetric Gaus- 
sians for signal and continuum background events, and 
with a Gaussian for the combinatorial BB background. 
For the peaking BB background we use a Gaussian dis- 
tribution for the DK mode. For the DK* mode, an 
asymmetric Gaussian is used for the K'K'k mode and a 
sum of two asymmetric Gaussians for the i^gTr mode. 
The shape parameters of the threshold function for con- 
tinuum background are determined from data. All other 
PDF parameters are derived from the simulated events. 
In the fits we fix the numbers of peaking BB back- 
ground events, which are estimated from the PDG 
branching fractions [llj and MG efficiency evaluations. 

The number of signal events determined by the fit 
(A^sig) is used to calculate the branching fraction as 



6(5^ 



D^K 



On 



A^si 



A^B+ 



Bu ■ Bko 



where Nb+ is the total number of charged B mesons 
in the data sample (equal to the total number of all 
BB pairs produced, since we assume equal production 
of B+B- and B°B°), Bd and 6^0 are the branch- 
ing fraction for each D meson decay channel and for 
K^ — > TT+TT^ respectively [ll|, and Csig is the recon- 
struction efficiency for each D decay channel evaluated 
from MC events. The expression for B{B+ -> D+K*^) 
is obtained replacing Bj^o /2 with the branching fraction 
of K*'^ -> iC+TT", Bji^o. The likehhoods for individual 
channels are combined to derive average branching frac- 
tions for B+ -^ D+K° and B+ -^ D+K*°. 

The fit procedure is validated using an ensemble of 
simulated experiments (toy MC studies) with all yields 
generated according to Poisson distributions. The non- 
floating parameters of the fits as well as the shapes of 
the background threshold functions are fixed to the val- 
ues obtained from the MC samples. We define the pull 
for a variable x as the difference between the fitted xat 
and the mean generated value (a; gen), divided by the error 
CTcrr, a^puii = (xfit - (cCgon) )/o-orr ■ We use the negative er- 
rors for fitted values that are smaller than the generated 
ones and the positive errors in the opposite case. The 
procedure gives Gaussian-like pull distributions for each 
channel and thus no biases of the fit model were found. 
In Table IIIII we show resulting expectations of asymmet- 
ric errors for each channel. The 95% probability ranges 
for these errors obtained from toy MC studies are also 
shown. Tests of the fit procedure performed on the full 
MC samples give values for the yields compatible with 
the generated ones. 

The main results of the fit to the data are reported in 
Table HVl which gives the values of the fitted parameters 
for each D channel and for the combination of fits. The 
background yields are close to the expectations and the 



errors obtained on the branching fractions are in good 
agreement with the values reported in Table IIIII The 
leading contribution (as expected) is obtained from the 
Kirn mode. Likelihood fit projections of the toes and J^ 
distributions are shown in Fig. [2 In Fig. [3] we also show 
for illustrative purposes the fit projection for TOeSi after 
requiring J^ > 0, to visually enhance any possible signal. 



SYSTEMATIC ERRORS 

We consider various sources of systematic error. One 
of the largest contributions comes from the uncertain- 
ties on the PDF parameterizations. To evaluate the con- 
tributions related to the toes and T PDFs, we repeat 
the fit varying the MC-obtained PDF parameters within 
their statistical errors, taking into account correlations 
among the parameters (labeled as "PDF - MC" in the fi- 
nal list of systematic error sources) . Differences between 
the data and MC (labeled as "Data - MC PDF shapes" 
in the final list of systematic error sources) for the shapes 
of Toes and J^ distributions are studied for signal com- 
ponents using data control samples. i?° -^ D+tt^ and 
B'^ — > D'^p" selected events are used to obtain the jties 
and J^ parameters for the DK and DK* modes, respec- 
tively. The analysis strategy is the same as for the signal 
events except for specific criteria to select K'^ or K*'^. For 
the continuum background, we estimate this uncertainty 
by repeating the fit using the PDF parameters obtained 
from off-resonance data instead of those from continuum 
MC. Finally, for the BB background, we estimate this 
uncertainty by leaving the parameters that describe the 
BB combinatorial background as free variables in the fit 
(separately for toes and J^). The systematic uncertainty 
is defined as the difference in the branching fraction re- 
sults from the nominal and alternative fits summed in 
quadrature. 

The systematic errors on the signal reconstruction effi- 
ciency include the uncertainty due to limited MC statis- 
tics, uncertainties on possible differences between data 
and MC in tracking efficiency, A'^ and 7r° reconstruction, 
and charged-kaon identification. In addition, there are 
additional contributions to these uncertainties originat- 
ing from the disagreement between data and MC distri- 
butions for all the variables used in the selection. These 
are estimated by comparing the data and simulation per- 
formance in control samples. To evaluate the uncertain- 
ties arising from peaking background contributions, we 
repeat the fit by varying the numbers of these events 
within their statistical errors. The uncertainties on the 
branching fractions of the sub-decay modes are also taken 
into account. The uncertainty on Ngg (1.1%) has a neg- 
ligible effect on the total error. 

The systematic uncertainties on the branching frac- 
tions are summarized in Table [V] All the uncertainties 
are considered to be uncorrelated and are treated sepa- 
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FIG. 2: (color online) Projections of the 2D likelihood function onto the rriES (top two rows) and J- (bottom two rows) axes for 
(a) Kn-K, (b) Ktttvtt^, (c) KItt and (d) K'^tt-k° for the B+ -^ D+ K^ mode, and (e) K-k-k and (f) K^-k for the B+ -^ D+ K*° 
mode. The data are indicated with black dots and error bars and the (blue) solid curve is the projection of the fit. 




FIG. 3: (color online) From top left to bottom right: tties projection for (a) Ktttv, (b) KTVTVTr'\ (c) K^n, and (d) KgTnr'^ for 
the B^ — i- D^ Ks mode and (e) Ktvtt and (f) K^tv for the B'^ — > D^ K*^ mode. The data are indicated with black dots and 
error bars and the different fit components are shown: signal (black solid curve), combinatorial BB (green dotted), continuum 
(magenta dot-dashed) and BB peaking background (red dotted) and the blue solid curve is the projection of the fit. We require 
J^ > to visually enhance the signal component. Such a cut has an approximate efhciency of 70% for signal, while it rejects 
more than 80% of the continuum background. 
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ratcly for each channel. 



RESULTS FOR BRANCHING FRACTIONS 

The final likelihood for each decay mode is obtained 
by convolving the likelihoods for the measured branching 
fractions with Gaussian functions of width equal to the 
systematic uncertainty. 

The final results including systematic uncertainties are 

B{B+ -^ D+K^) = (-3.8tli) X 10-^ 
B{B+^D+K*^) = (-5.3±2.7) X 10"^ 

Since the measurements for the branching fractions 
are not statistically significant, following a Bayesian ap- 
proach and assuming a flat prior distribution for the 
branching fractions, we integrate over the positive por- 
tion of the likelihood function to obtain the following 
upper limits at 90% probability: 

B{B+ -^ D+K") < 2.9 X 10~^, 
B{B+ -^ D+K*") < 3.0 X 10"^ 

The -B"*" — >■ D^ K'^ result represents an improvement 
over, and supersedes, our previous result [7[, while the 
5+ -^ D^ K*^ result is the first for this channel. 



CONCLUSIONS 

In summary, wc have presented a search for the rare 
decays B^ -^ D^K^ and i?+ -^ D^K*^ ^ which are pre- 
dicted to proceed through annihilation or rescattering 
amplitudes. We do not observe any significant signal and 
we set 90% probability upper limits on their branching 
fractions. 
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TABLE V: Systematic errors on branching fractions for B'^ — >■ D^ K'^ and B^ — >■ D^ K*'^ decay channels. All quantities are 
given in units of 10~ . 







B+ - 


-+ D+K'' 


B+ ^ D+K'° 




K-K-R 


KTVTnv" 


K°s7V Kl7vn° 


K-KTT K°sn 


PDF - MC 


+0.8 
-0.8 


+6.2 
-3.4 


+5.3 +7.3 
-4.4 -8.8 


+0.6 +3.1 
-0.9 -3.6 


Data-MC PDF shapes: 



Continuum background 0.2 0.4 1.4 0.5 0.1 1.7 

SB background 0.7 1.6 2.5 5.0 1.0 4.4 

Signal < 0.05 9.2 5.6 0.9 0.9 3.1 



Efficiency error: 


Reconstruction efficiency (MC) 0.1 
Data-MC 0.2 


0.6 

0.8 


<0.05 
<0.05 


0.9 
0.5 


0.1 
0.2 


0.5 
0.3 


Peaking background < 0.05 
B errors 0.3 


0.5 
0.3 


0.2 
<0.05 


0.2 
0.4 


< 0.05 

< 0.05 


0.1 
0.1 


Combined tH 


+ 11.3 
-11.8 


+ 8.2 
-9.3 


+9.0 
-12.5 


+ 1.5 
-1.8 


+6.4 
-7.4 



